Error threshold estimates for surface code with loss of qubits 
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We estimate error thresholds for surface code in the presence of loss via an analytical method 
developed in statistical physics. The error threshold for the surface code is closely related to a special 
critical point in a finite-dimensional spin glass, which is disordered magnetic material. We compare 
our estimations to the heuristic numerical results reported in earlier studies. Further application 
of our method to the depolarizing channel, a natural generalization of the noise model, unveils its 
wider robustness even with loss of qubits. 



Introduction — Against corruption by environmental 
noise as well as imperfection in implementation, the state 
of qubits describing quantum information cannot be sta- 
ble and must be recovered by elaborated procedures, 
quantum error correction [J f Quantum error correc- 
tions usually work on the computational error on qubits, 
which do not go out of the basis for computations. There- 
fore errors come from losses of the physical resource, 
qubits, can deprive of the performance of error correc- 
tion. However, if one can detect and identify locations of 
losses, a modified scheme can recover the original infor- 
mation. Stace and Barret have suggested a topological 
error-correcting code, a family of Kitaev's surface codes 
■3j, which is robust against both of the computational 
errors and losses by a modified scheme according to the 
location of the lost qubits Q • 

In the present paper, we estimate precise values of the 
error thresholds for the modified error correcting code 
against both of the computational errors and loss of 
qubits by use of a systematic theory developed in sta- 
tistical physics. The key of our analysis is hidden in the 
disordered magnetic system, spin glasses. Several spin 
glass models have a special symmetry with exact solv- 
able subspace known as Nishimori line 0, Q . The critical 
point in this subspace, termed as the multicritical point, 
corresponds to the optimal error threshold in the topo- 
logical quantum error correcting code Q • A combination 
of the duality with the real-space renormalization tech- 
nique, which are often used to identify the singular points 
in statistical mechanical models, can derive the precise 
estimations for the optimal error thresholds 
systematically approach the exact solutions 
use of this method, 
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we fill the blank on the analytical 
study for the optimal error thresholds on several topolog- 
ical error correcting codes with loss of qubits. The results 
reported in this paper provide upper bounds against er- 
ror rates for any error-correcting schemes. They serve as 
important benchmarks with which any constructive error 
correcting procedure as recently proposed in Ref. [7| can 
be compared. 

Surface code and spin glass — Let us consider qubits 
set on each edge (ij) of the square lattice embedded 



on a torus (genus 1). We define the star operator 
Xuj\ for each site s, and plaquette operator 



Zuj\ for each plaquette p (site on the dual 




lattice), where X and Z are Pauli matrices. The product 
consists of four edges adjacent to each site or plaquette. 
The stabilizer group is given by the simultaneous eigen- 
states with the positive eigenvalues for these operators 
X s and Z p . Since the star and plaquette operators con- 
sist of unit loops on the dual and original square lattices, 
any contractible loop by X s and Z p products on each lat- 
tice acts trivially on the codespace. On the other hand, 
any non-contractible loops on the lattice can map the 
codespace to itself in a nontrivial manner. If we set Lx L 
lattice on a torus, we have 2L 2 qubits and 2(L 2 — 1) stabi- 
lizers. The remaining degrees of freedom of 2 implies ex- 
istence of two non-contractible loops, winding around the 
hole of the torus L v and winding around the body of the 
torus L t , and ones L* and L* t on the dual lattice. These 
loops can be written in terms of the products of operators 

as Z v = Il(y)ei v Z (v)> Xv = H{i])eL* v X (ij)> Z t' and X u 
which are termed as logical operators. The logical oper- 
ators can form Pauli algebra of two effective qubits en- 
coded in the topological degrees of freedom on the torus 
as [Z v ,Z t ] = [X v ,X t ] = 0, and X t Z v = {-\) s ^Z v X t . 
The combinations of non-contractible loops yield 2 4 = 
16 different homology classes for the original and dual 
square lattices on a single torus. We need to distinguish 
them for protecting the information from corruption. 

In order to evaluate the performance of the error- 
correcting code, let us define a noise model where each 
qubit independently gets errors as 

P^Pip+ (pxXpX + PY YpY +p z Z P Z). (1) 

Although, if we employ the following analytical method, 
we can estimate precise values of the error thresholds for 
"any" cases of pi,px ,Py , and pz, we restrict ourselves 
to two cases: px = Pz — P, Py — p 2 and p/ = (1 — p) 2 
(uncorrected case), and px — Py — Pz = p/3 and 
pi = 1 — p (depolarizing channel case) for simplicity, 
where < p < 1. The error Yuj\ can be regarded as a 
multiple error Xuj-\ and Zuj\ . The errors Zuj\ and Xuj\ 
can be described as chains E and E* on the original and 
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dual lattices. The endpoints of the error chains dE and 
dE* can be detected by applications of star and plaque- 
tte operators due to anti-commutation of adjacent errors 
with operators. From the knowledge of endpoints dE 
and dE* without the homology class of the error chains, 
error syndrome, we infer the most likely homology class 
of error chains, while considering any reasonable choices. 
Since E' = E + C and E'* = E* + C*, where C and C* 
are the contractible loops on both of the lattices, are in 
equivalent class with the error chains, the probability for 
the homology class E and E* of the error chains can be 
written as [3j 

P(E,E*\dE,dE*) = P(E,E*)/^2,P Di {E,E*), (2) 

i 

where P(E,E*) cx £ cc , JJm exp^rg + Kt?' rf/) 
for the uncorrelated case. The summation is taken over 
all the possibilities of C and C* , and the product is 
over all the edges. The parameter K stands for the im- 
portance/preference to choose the inferred error chain. 
The quantity in the denominator Pn i (E, E*) denotes the 
probability with the different homology class specified by 
the logical operators Di (i = 1, 2, • ■ ■ , 2 4 ). We here use 
to represent the inferred error chains, which takes ±1 
(t? < 0, when (ij) £ E), and also for E* , C and C* . 

The loop constraints fJto) T f 3 = 1 and Yluj) T ij* = 1 al- 
low us to use another expression by the Ising variables 

T W) = ai<T j' an< ^ T (ij) ~ a i a j ^ or eacn lattice on the 
torus. By use of these expressions, we can find that 
P(E,E*) is written as square of the partition function 
of the ± J Ising model 

P(E,E*) cx ^2 Y[e K ( T 5 a ^ +T v a >'j\ (3) 

<J,a" (ij) 

where and t~ are the signs of the quenched ran- 
dom couplings in context of spin glasses. When we set 
K = Ki n d., where exp(2if inc i.) = (1 — p)/p (Nishimori 
line), the inference of the error chains is an optimal re- 
covery procedure to identify the most likely homology 
class |5j . Each of the quenched random couplings follows 
the distribution function of the error chains P(E, E*) = 
Tltij) P( T ^j)P( T ij ) f° r the uncorrelated case, where 

P(r|) - (l-p)8 Til>1 +pS Til> _ 1 . (4) 

Similarly, we can evaluate the probability P(E,E*) for 
the homology class of the error chains E and E* for the 
depolarizing channel case as 

P(E,E*)*]Je KT $+ KT v +KT 3 T %' , (5) 

(ij) 

where we set the parameter K = -Kdep. as on the Nishi- 
mori line exp(4iidcp.) = 3(1 —p)/p. This is written in 



terms of the partition function of the eight- vertex model 
with quenched random interaction Q 

Zdcp . = j2Y[e K ^^+^' ->;+^-5' ( 6 ) 

a, a* {ij) 

where tE and Ty follow the distribution function 
through P(E,E*) = U {ij) P d ep.(r| ,-rff) asP dep .(l, 1) = 
l-p while P dep .(l, -1) = Pde P .(-l, 1) = Pde P .(-l, -1) = 
p/i. We emphasize that, if we tune the probability func- 
tion appropriately, we can apply our analysis as shown 
below to inhomogenius case with px ^ Py ^ Pz- 

In context of the statistical physics, Di represents the 
domain wall. In the low-temperature region implying 
a small p, the order of the degrees of freedom sup- 
presses the fluctuation of the domain wall. The cost 
for free energy difference due to the domain wall di- 
verges as J2 E ,E* P ( E > E *) P (E'E*\ dE > dE *) -> 1 for 
L — > oo. This means that we can infer the equiva- 
lent class with the original error chains. On the other 
hand, in the high-temperature region, the cost vanishes 
and J2e.e* P{E,E*)P{E,E*\dE,dE*) -> 1/16. This 
implies that the failure of the recovery occurs at the crit- 
ical point. Therefore the location of the critical point on 
the Nishimori line, the multicritical point, identifies the 
optimal error threshold. 

Loss of qubits and bond dilution — Loss of qubits on 
the lattice implies the modification of the stabilizers as 
well as the logical operators. However we can reform a 
complete set of stabilizers even on the damaged lattice 
due to loss of qubits following the proposed scheme in 
Ref. [H . The effect of lost qubits appears in the pattern 
of error chains E and E* , and their weight for the proba- 
bility, which degrades the performance of error-correcting 
code. To infer the most likely homology class based on 
the knowledge of the error chains on the damaged lattice, 
we reconstruct the original lattice by assigning of weight- 
zero edges on the lost qubits and irregular weight edges p' 
adjacent to the lost qubits as 1 — 2p' n , where n is the num- 
ber of the shared qubits in adjacent edges as in Fig. [T] 
The weight-zero edges imply that we need to consider a 
diluted version of the original spin glass system as in Eqs. 
d2]) and It can be achieved by a simple modification 
of the distribution function for Ty and into, for the 
uncorrelated case, P q (r^) = (1 — q)P(r^) + q8 T E and 
P q (r^ ), where q denotes the ratio of loss of qubits. Sim- 
ilarly, for the depolarizing channel case, P^ ep (r^ , Ty ) = 

(1 - q)P dep .( T ij^ T ij ) + IKst?;* ,0' In addition, we have 
to take into account effects of irregular weight edges p' 
adjacent to the lost qubits as carefully discussed in Ref. 
9] . The effect can be described by highly correlated dis- 
tribution function depending on the pattern of the lost 
qubits, although we omit its detailed expression. 

Duality analysis for spin glasses — Analyses to clarify 
the critical phenomena in finite-dimensional spin glasses 
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FIG. 1: Reconstruction of the damaged lattice by use of the 
weight-zero and irregular weight. The dashed line denotes the 
lack of qubits and weight-zero edge. The bold line expresses 
the irregular weight edges after the reconstruction. In this 
case, n = 2. 



are intractable in general. However a recent develop- 
ment in the spin glass theory enables us to estimate the 
precise value of the special critical point on the Nishi- 
mori line, which corresponds to the optimal error thresh- 
old [10Ml3j . The method as shown below is based on 
the duality, which can identify the location of the criti- 
cal point especially on two-dimensional spin systems [3] ■ 
Let us review the simple pure Ising model case at first. 
The duality is a symmetry argument by considering the 
low and high-temperature expansions of the partition 
function Z = Yiuj) exp(ifo"j<7j). The painful cal- 
culation of both expansions can be replaced by a sim- 
ple manipulation with the binary Fourier transformation 
for the local part of the Boltzmann factor, namely edge 
Boltzmann factor xo — exp(K) and x\ = exp(— K) [14j . 
The low-temperature expansion can be expressed by xq 
and X\. On the other hand, the high-temperature ex- 
pansion is given by the binary Fourier transformation 
Xq = (xo + xi)/V2 and x* = (xq — Xx)/y2. We use this 
fact and find a double expression of the partition function 
as 

{x,{K)Y L2 z( Ul {K)) = {x*(K)r L2 z(ul(K)), (7) 

where z is the normalized partition function z{u\) — 
Z/{x Q (K)} 2L2 and z(u\) = Z / {x* q (K)} 2l2 . We here de- 
fine m{K) = xi(K)/x (K) = exp(-2if) and = 
x*{K)/xq{K) — tanhif. The well-known duality rela- 
tion exp(— 2K*~) = tanhK is given by rewriting u*(K) by 
u\{K*), which implies a transformation of the tempera- 
ture. Then the principal Boltzmann factors xq{K) and 
Xq(K) with edge spins parallel holds xq(K c ) = Xq(K c ) at 
the critical point exp(— 2K C ) — t&nhK c . 

We employ the replica method, which is often used in 
theoretical studies on spin glasses, in order to generalize 
the duality analysis to spin glasses 1(J 11|. Let us con- 
sider the duality for the replicated partition function as 
[Z£ d ] and [ZdepJ simply [Z n ], where [• • •] is the configu- 



FIG. 2: Clusters for the uncorrelated case (the left three pan- 
els), and depolarizing channel case (the right two panels). The 
cluster Boltzmann factor is defined without omitting many 
body interactions generated after the decimation of the black 
spins on the original lattice. The dotted line denotes the dual 
lattice for a* and tE^. The subscript k denotes the configu- 
ration of the white spins. 



the replicated partition function as 

{x (q 1 K)} 2L2 z(u 1 (q,K),u 2 (q,K),--- 
= {x* (q,K)} 2L2 z(u* 1 (q,K),u* 2 (q,K)^ 



(8) 



where the subscript of and u\ stands for the number 
of anti-parallel pair among n replicas on each edge. Un- 
fortunately we cannot replace u* k (q,K) by Uk(q*,K*) as 
the pure case, since the replicated partition function is 
multivariable. Nevertheless we can estimate the precise 
location of the critical point even for spin glasses by con- 
sidering a wider range of the local part of the Boltzmann 
factor given after the summation of the internal spins. 
For instance, in the case on the square lattice, we define 

(s) 

the cluster Boltzmann factor x k , where the subscript 
k denotes the configuration of the edge (white-colored) 
spins and s expresses the size of the cluster in Fig. [2j We 
set the equation to lead the location of the critical point 
as, inspired by the case without quenched randomness, 
x (K)=x* (K) EMI, 



4 s) (q,K)=x^'(q,K) 



0), 



(9) 



rational average for the quenched randomness according 
to the distribution functions. The multiple (2 n ) Fourier 
transformation again leads us to the double expression of 



The equality even for s — (edge) can give the precise 
solutions of the critical point for the multicritical point of 
±J Ising model q = as p c = 0.1100 HEED]- Although 
the above method is not exact, if we increase the size 
of the used cluster, we can systematically approach the 
exact solution for the critical points of the ± J Ising model 
in the higher temperature region than the Nishimori line 

Results — In the present study, we consider two sizes 
of the clusters as s — 1 and 2 as well as a single edge 
(s = 0) for the uncorrelated case and a single cluster 
(s = 1) and a single crossing edges (s = 0) for the depo- 
larizing case as in Fig. [2j We show several estimations 
given by Eq. for the uncorrelated case in Table. Q] 
Although, for the uncorrelated case by s = 2 cluster, we 
have considered the highly correlated distribution func- 
tion by taking into account the effect of the irregular 
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weight, the results have not been changed from those by 
a simple distribution function P q (E), which is the same 
as one for the bond-diluted spin glass. As discussed in 
Ref. [9( , the highly correlation between the loss of qubits 
is found to emerge as a finite-size effect in the numerical 
investigation (for q > 0.45). Such the complicated effect 
does not spoil our analysis. All the results for any q does 
not show drastic changes dependently on the size of the 
used cluster. It means that our analyses are enough cor- 
rect to estimate the exact locations of the optimal error 
threshold. The optimal error thresholds indicate the up- 
per bounds for error threshold by any heuristic methods. 
As shown in Fig. [3j we compare our results with the 
inference by use of the matching algorithm namely, the 
ground state as K — > oo [4j , in which we denote the error 
thresholds as p° c . We confirm that the heuristic matching 
algorithm of inference gives p c « p° c , presumably p c = p a c 
for large q. We also give several results for the depo- 
larizing channel case in Table. HU Similarly to the case 
without loss of qubits {q = 0) as have been reported in 
Ref. [8(, the depolarizing channel is more resilient than 
the uncorrelated case even with loss of qubits. For com- 
parison, let us take an earlier study on an error recovery 
procedure for the depolarizing channel in Ref. Q. Our 
result implies that there still possibility to improve the 
performance of such a constructive procedure. All the 

q p c = 0) p c (s= 1) p c (s = 2) p° c H] 



0.00 


0.11003 


0.10928 


0.10918 


0.10486 


0.10 


0.09240 


0.09196 


0.09189 


0.08816 


0.20 


0.07245 


0.07235 


0.07233 


0.06997 


0.30 


0.04984 


0.05004 


0.05009 


0.04836 


0.40 


0.02462 


0.02492 


0.02500 


0.02561 


0.45 


0.01155 


0.01174 


0.01179 


0.00757 




TABLE I: Comparison of the approximations for s = 0, 1, 
and 2 for the uncorrelated case and by a heuristic method [J . 



FIG. 3: (Color online) Results for the uncorrelated case. The 
dots stand for numerical data given in Ref. Q, where the 
error bars are suppressed in this scale. The curves almost 
overlapped in this scale express our results by the duality. 



loss of qubits, via a finite-dimensional spin glass the- 
ory, for both of the uncorrelated and depolarizing chan- 
nel cases, and shown more resilience of the depolarizing 
channel even with loss of qubits. 

In the sense of study on spin glass, the comparison 
between the error thresholds p° c by a suboptimal method 
corresponding to the inference in the ground state [4j and 
optimal ones p c shows a fascinating feature of the phase 
boundary of the ±J Ising model as p c ~ p® 5]. The 
future study will be desired for solving the remaining 
problem on a realm of spin glasses: p c = p° c or not. 
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q 


Pc (s = 0) 


p c (s = 1) p c \7} 


0.00 


0.18929 


0.18886 0.164 


0.10 


0.16025 


0.15985 


0.20 


0.12690 


0.12656 


0.30 


0.08844 


0.08819 


0.40 


0.04454 


0.04440 


0.45 


0.02121 


0.02114 



TABLE II: Results for the depolarizing channel case. 

obtained values are almost stable in the third digits. In a 
practical sense, our estimations for error thresholds serve 
as the reference values. 

Conclusion — We have estimated the error thresholds 
for the topological quantum error-correcting code with 
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